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SECTION  1 
INTRODUCTION 


Under  the  sponsorship  of  the  Defense  Nuclear  Agency,  a  quantitative  methodology  was 
developed  to  treat  identified  ermrs  affecting  targedng  predictions.  Described  herein  is  a  derivative 
of  that  methodology  specifically  developed  to  support  the  SEPW  Phase  2  Study. 

This  targeting  uncertain  ies  methodology  quantifies  the  effect  of  systematic  errors  on 
damage  expectancy,  DE  (a  measure  of  effectiveness  used  by  the  JSTPS  to  quantify  the  probability 
of  achieving  at  least  the  specified  level  of  target  damage).  Tlie  methodology  evaluates  and  ranks 
the  effect  of  uncertainties  for  applications  involving  a  single  buried  target  and  a  single  earth 
penetrating  weapon.  The  methcxlology  models  the  penerration  of  the  EPW,  the  ground  shock 
environment,  and  the  probability  of  damage  to  the  underground  target.  Potential  error  sources 
treated  in  the  methodology  are  impact  velocity,  geology,  depth  of  burst,  target  hardness,  ground 
shock  prediction  model,  ground  shock  depth-to-effect,  weapon  aiming  accuracy,  distance  damage 
prediction  model,  and  target  depth,  location  and  size. 

The  methodology  makes  minimal  use  of  assumed  distributions  and  Monte  Carlo 
techniques;  instead  it  uses  three  points  to  describe  the  inherent  dispersion  for  all  but  one  error 
source.  DE  is  calculated  for  each  combination  of  error  source  values.  The  values  are  sored  and 
used  to  define  a  distribution  for  DE.  In  addition,  an  analysis  of  variance  technique  is  used  to 
calculate  the  relative  contributions  of  each  error  source  to  the  uncertainty  in  DE.  The  added  benefit 
of  ranking  the  error  sources  in  terms  of  their  contribution  to  the  uncertainty  in  DE  can  provide 
insight  into  deeply  buried  targeting  problems. 

The  uncertainties  methodology  is  PC-compatible,  fast  running,  and  provides  results 
which  compare  favorably  with  more  rigorous  mathematical  methods. 

The  methodology  will  be  described  in  Section  2,  followed  by  an  example  problem  in 

A 

Section  3.  Section  4  provides  a  summary.  Hardware  requirements  and  runtime  procedures  are 
given  in  Section  .5.  Section  6  provides  references. 
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SECTION  2 

UNCERTAINTIES  METHODOLOGY 


The  uncertainties  methodology  for  deeply  buried  targets  is  a  PC-based  algorithm  which, 
unlike  a  Monte  Carlo  based  routine,  makes  minimal  use  of  assumed  distributiom\.  Instead  it  uses 
three  points  to  describe  the  inherent  dispersion  of  all  but  one  of  the  identified  error  sources.  This 
allows  infinite  flexibility  in  describing  the  suspected  error.  The  three  points  may  be  thought  of  as  a 
nominal  value  and  lower  and  upper  bounds;  mean  plus  or  minus  one,  two,  or  three  standard 
deviations;  median  and  quartiles;  or  for  best  agreement  with  Monte  Carlo  methods,  5th,  50th,  and 
95th  percentiles.  Generally,  all  combinations  of  values  for  the  error  sources  are  propagated 
through  the  methodology,  creating  a  distribution  of  outcomes  and  allowing  for  an  analysis  of 
variance  (ANOVA)  technique  to  calculate  the  relative  contribution  of  each  error  source  to  the 
uncertainty  in  outcome.  A  more  complete  di.scussion  of  the  3-point  concept  and  its  use  is  given  in 
two  references  (Ref.  1)  and  (Ref.  4), 

Potential  error  sources  treated  in  the  methodology  include  EPW  impact  velocity,  geology, 
depth  cf  burst  (DOB),  target  hardness,  ground  shock  prediction  model,  ground  shock  depth-to- 
effect,  target  depth,  target  location,  target  size,  weapon  aiming  accuracy,  and  distance  damage 
prediction  model.  An  overall  view  of  the  uncertainties  methodology  is  given  in  Figure  1 . 

As  shown,  the  uncertainties  methodology  is  performed  in  two  par's.  The  first  part 
(DTEUA)  analyzes  the  uncertainty  in  depth-to-effect  (DTE),  given  uncertainties  in  EPW  impact 
velocity,  target  site  geology,  target  hardness,  and  the  DUG  1C  ground  shock  mode,'.  Outputs  to  the 
DTE  analysis  include  uncertainty  in  DTE,  and  relative  contributions  of  each  errer  source  to  the 
uncertainty  in  DTE.  In  addition,  a  value  for  PSP  (probability  of  successful  penetration)  is 
calculated  based  on  the  number  of  Monte  Carlo  trials  which  exceeded  recommended  levels  of  axial 
loading,  or  in  which  the  penetration  was  sufficiently  shallow  to  incur  a  rebounding  condition. 

The  second  part  of  the  uncertainties  methodology  is  concerned  with  uncertainties  which 
influence  DE  (and  were  not  examined  in  the  first  part  of  the  methodology).  Error  sources  of 
interest  include  target  size,  target  depth,  CEP  (weapon  aiming  accuracy),  target  location,  slant 
range  damage  sigma,  and  DTE  (which  was  calculated  in  the  first  pan  of  the  methodology).  Values 
for  these  error  sources  are  used  as  inputs  to  the  DUGDE  code  along  with  the  PSP  value  calculated 
in  DTEUA.  Outputs  from  the  DE  analysis  include  a  best  estimate  for  DE  as  well  as  one  standard 
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deviation-type  uncertainty  bounds  for  DE.  Also  calculated  are  the  relative  contributions  for  each  of 
the  error  sources  to  the  uncertainty  in  DE. 
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Figure  1.  Overview  of  the  uncertainties  methodology  for  deeply  buried  targets. 


2.1  DTE  ANALYSIS. 


This  part  of  the  methodology  is  concerned  with  the  penetration  of  the  EPW  into  the  target 
site  geology,  propagation  of  the  ground  shock  down  toward  the  target,  and  the  calculation  of  the 
depth- to-effect.  A  logic  diagram  of  the  DTE  analysis  is  shown  in  Figure  2. 
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Itu’  I M  I  rthrtlvii^  i'xrtnunoN  thr  onci  i  o|  iiiu'i*i!rtiniv  tti  itnpiu'i  v«'l«H'liy,  jjc*olouy,  tiwyi'i 
ItiihlnoM,  iMui  ot«  t)!)'  Ihirc  |Htinu  lur  tiHcd  lo  ;h('  iIInpoimod  lit 

imp.u » iiv,  (rtiufl  m  omu  •>!  !»<*rtk  luut  v  oiilulriu^?  ii\  iiroumi  slun  k 

MUsU't  lOSklltS 

Hosst^vci  It)  uiiliM  »o  assess  iho  unpas  i  of  iiiuenaiiuy  In  )ienli»ay  on  tleplh  of  penoiraiton 
aiul  thus  DIP,  a  Moiue  ('aii»>  levhnupio  is  useil,  Uiuenainty  is  as.sljpu'il  loeiich  layer's  iliUkness 
aiul  pi'iieiialnltiy  i  liven  an  iinpaei  veUn  Hy.  a  Ntoiue  i'aiK»  sampling  i>*i  ht",ipie  is  usetl  lo  sample 
iMvIi  lavi'i  s  iliukness  iassniue<l  lo  he  uniformly  illsfrilupetll  ami  peneirahilliy  (assnmtii  lo  l*e 
iioimallv  itisiiihnieil)  Tlien  ilie  peneiration  equations  are  used  to  ealciilate  the  tiepili  of  peneinuion 
inio  Iho  sampled  jjeolo^y  ami  the  maximum  deceleration  es|X‘nem'ed  hv  ihe  peneiriiior, 

Hiese  ivvo  steps  (sample  jieologs,  cahulaie  peneinuion)  are  repeainl  for  a  yiven  numkM 
o.|  uials  in  Ollier  lo  achieve  a  distiihuiion  of  penetration  dcptlv*.  'vhich  lepresents  the  unceriainty  in 
Dt  )H  .lue  lo  ihe  uiu  eriamiy  in  ueolojjy,  Prom  this  distrihuiion,  the  5ih.  SOih,  and  ‘>.'Sih  percentiles 
■le  read  and  saved  tis  inputs  for  the  ground  sluK'k  nuxiel.  In  addition,  the  values  calculatcil  for 
depth  of  penoiiaiion  and  maximum  deceleration  for  each  trial  are  checked  against  relxuimling  ami 
axial  loading  ci  iiena  to  determine  the  piuhahiliiy  of  successful  peneinuion  (PSP).  The  value  for 
PSP  IS  usoil  as  input  to  the  DP  analysis, 

['he  caicul  uioii  of  uncertainly  in  DOH  due  to  uncertainly  in  geology  is  repeated  for  each 
impact  vchK’iiy,  Thus  if  Ixuh  geology  and  impact  velocity  arc  assigned  uncertainties,  nine  DOBs 
will  he  calculated. 

I'acli  DOB  is  used  as  an  input  to  the  ground  slux'k  model,  which  predicts  the  down-axis 
depth  versus  peak  stress  profiles.  Values  for  depth  to  each  target  iuirdness  are  detennined,  us  well 
.IS  the  uncertainly  in  those  depths  due  to  potential  rmxieling  errors.  Since  there  lu^  nine  DOBs  and 
nine  depths  (I.)TP.s)  read  from  the  ground  shixk  model  output  for  each  DOB,  in  all  81  (3‘^)  DTEs 
arc  calculatcil.  i'hcse  81  D  TE  values  ore  sorted  to  form  a  cumulative  density  function  (cdO  for 
Dl'E, 


5 


l'olU>wtni£  illxcussion  of  ihe  pchetniiion  mut  shiKk  tiunJels,  ihc  c«lculntlt>n  of  tin* 

n*l(»iivi'  i'ontiil)»tit)Hs  of  rtn  h  of  the  four  crr\)r  stnirtYs  will  he  tliscusscd,  «,<  well  nit  the  coiwersloit 
of  the  K1  cull  ulitieil  l)  ri\s  ittio  the  1  point  whente  for  input  into  the  DU  tuuklysi,*), 

^,1,1  IVncinuion  MtHlel, 

The  tnotlel  used  to  detenninc  the  depth  of  penoivution  of  the  UPW  is  bused  on  the 
eiiuiuit>ns,  definitions,  und  ulgorlthms  in  two  tlcvclopincnt  reports  by  C.W,  Young  (Ref.  b,  Ref, 
7),  I'he  nttxlel.  culled  PHNDUPTII,  is  documented  in  its  user's  munuul  (Ref.  5).  The  nuxlel 
calciilutes  the  depth  of  penetration  of  un  UPW  through  u  layered  geology.  Clencrully,  the 
pcnctrutor's  vekwity  e.xitin}i  a  layer  is  calculated,  then  applied  as  an  impact  vckx'ity  into  the  next 
layer,  litis  is  repeated  for  each  layer  until  the  |)enetraior  stops. 

Inputs  to  the  tmxlel  include  those  associated  with  both  the  weapon  und  the  geology. 
Weapon  inputs  ittcludc  impact  angle  und  velocity;  pcnetnitor  weight  and  diameter,  nose  shape,  nose 
length,  and  perfommnee  coefficient.  Geology  inputs,  for  each  layer,  include  a  layer  type  (i.c,. 
soil,  rock,  concrete,  or  froren  soil),  a  thickness,  and  a  penetrability  number  (called  the  S-number), 
The  harder  a  layer  is.  the  lower  its  S-number  will  be,  and  in  general  the  shallower  the  penetration 
depth, 

PUNDEl^H's  output  includes  the  path  length  and  vertical  component  of  the  depth  of 
penetration,  as  well  as  the  layer  in  which  the  penetrator  comes  to  rest.  Additional  outputs  include 
the  peak  and  average  axial  deceleration  experienced  by  tlie  penetrator  in  each  layer,  and  the  velocity 
as  the  penetrator  enters  each  layer. 

2.1.2  Ground  Shock  Model. 

The  values  calculated  for  depth  of  penetration  due  to  uncertainty  in  both  impact  velocity 
and  geology  are  used  as  input  to  the  ground  shock  model  along  with  penetrator  yield  and  geology 
information  (e.g.,  layer  thickness,  air-filled  voids,  density,  seismic  and  loading  velocities,  and 
attenuation  exponent).  Tlie  ground  shock  model  is  a  DNA-approved  engineering  algorithm,  called 
DUG  1C,  which  predicts  the  down-axis  (ID)  depth  versus  stress  profiles  beneath  near  surface 
explosions.  The  DUG  1C  computer  code  is  described  in  detail  in  its  user's  guide  (Ref.  2)  and 
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(«?chnic»il  inAnuH.1  ^Rcf,  In  uddition  to  down  axis  depths^iO'Crtccf,  DUO  1C  spccilles  «  geoi«eu7 
lor  consiniciing  iso-siitiss  conluurs,  us  shown  in  Figure  .V 

'n»e  DUGiC  ''odc  leproduncs  most  observed  contained  burst  explosion  duta  in  media 
rungittg  from  Ux)sc  soil  to  hard  rock,  and  is  accurate  to  a  fretor  of  two  in  stress  at  a  given  dcptlt. 


Ground 


Figure  3.  DUG  IC-dcfined  iso-stress  contour. 


Internal  to  DUGIC  is  a  coupling  curve  (see  Figure  4)  which  defines  the  portion  of  energy 
that  is  coupled  into  the  ground.  All  coupliiig  values  are  relative  to  1.0  for  a  fully  coupled  burst. 
Depending  on  where  the  DOB  is  on  the  coupling  curve,  even  a  small  uncertainty  in  DOB  can  have 
a  large  effect  on  the  coupling  factor,  which  is  directly  related  to  the  determined  DTE. 
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l.O 


Figiupc  4.  Coupling  curve  inherent  in  DUG  1C. 


Figure  5  integrates  the  penetration  and  ground  shock  models  to  show  tlie  effect  on  DTE 
of  aTi  uncertainty  in  DOB  (10  m  ±  a  factor  of  2)  for  a  1  MT  EPW. 


Figure  5.  Effect  of  DOB  uncenainty  on  DTE. 

* 
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2,1.3  Calculation  of  Relative  Contributions. 

Recall  that  81  (3“*)  DTEs  are  calculated,  where  each  DTL  value  is  from  a  different 
combination  of  values  for  each  error  source.  By  calculating  all  combinations  and  storing  which 
combination  of  inputs  ic.sulted  in  each  DTH  value,  an  analysis  of  vaiiance  (ANOVA)  technique  is 
used  to  determine  the  relative  contribution  of  each  error  source  to  the  uncertainty  in  DTE. 

The  ANOVA  methodology  is  comparable  to  a  four-way  ANOVA.  Two  i modifications  to 
tire  classical  ANOVA  technique  were  made  to  suppress  the  impact  of  potential  outliei  .  hich  might 
arise  through  use  of  only  three  points  to  describe  the  uncertainty  in  each  error  soui^?.  The 
methodology  uses  median  values  rather  than  means,  and  uses  sums  of  absolute  differences  nrher 
than  the  sums  of  squares.  The  methodology,  described  below,  was  originally  documented  c  • 
Binninger  and  Wright  (Ref.  1,  p  111-115). 

Tables  are  created  for  each  error  source  in  which  the  dispersion  in  each  aiw  is  due  only  to 
the  dispersion  in  that  particular  error  source.  Figure  6  shows  the  table  created  for  target  hardness. 


Row 

Geology 

Impact 

Velocity 

DUG 

Model 

DTE  when 

Target  Hardness  is: 

N  L  U 

1 

N 

N 

N 

-- 

. 

2 

L 

N 

N 

— 

-  - 

3 

U 

N 

N 

-  ■■ 

25 

N 

U 

U 

• 

• 

26 

L 

u 

U 

— 

-  - 

27 

U 

u 

u 

— 

—  — 

N:  Nominal  L:  Lower  U;  Upper 


Figure  6.  Table  used  to  calculate  relative  contribution  for  target  hardness. 


Contiibutions  from  each  row  are  summed  for  each  table,  e.g.,  for  target  hardness 


Sumn  ~  I  D'TEfow.N  ■  DTErow,L|  + 1  DTErow.u  ■  D  I'fcitow.N  | 

all  rows 


(1) 


Then  the  relative  contributions  are  calculated 


%H  = 


_ SumH _ 

{Sumiv  +  Sumo  +  Sumo  +  Sumn) 


(2) 


where  %h  is  the  percentage  ot  the  error  in  DTE  due  to  target  hardness,  Sumrv  is  the  sum  calculated 
for  impact  velocity,  f  ume  is  for  geology.  Sumo  is  for  DUGIC  modeling  error,  and  SumH  is  for 
target  hardness. 

When  examining  the  calculated  relative  contributions,  it  is  also  important  to  look  at  the 
corresponding  DE  uncertainty  bounds  (discussed  below).  For  example,  a  large  contributor  when 
the  DE  uncertainty  bounds  are  small  may  not  be  an  important  factor  if  its  contribution  decreases  as 
the  uncertainty  bounds  change.  Also,  relative  contributions  of  error  sources  whose  uncertainties 
may  be  decreased  by  research  should  be  weighed  against  contributions  from  those  error  sources 
whose  uncertainty  bounds  cannot  be  reduced. 

2.1.4  Determin.ation  of  Three  Points  to  Describe  DTE  Uncertainty. 

While  81  DTE  values  are  calculated  in  the  DTE  analysis,  the  next  part  of  the  analysis  (DE 
analysis)  requires  that  the  uncertainty  for  each  error  source  be  entered  in  the  3-point  format.  To 
this  end.  the  calculated  distribution  for  DTE  is  fit  to  a  log  normal  distribution  from  which  the  5th, 
50th,  ?nd  95th  percentiles  are  calculated  (Ref.  1,  p  134-140). 
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The  calculated  DTE  distribution  may  not  follow  any  known  distribution,  and  while  the 
5th,  50th,  and  95th  percentiles  could  be  read  from  the  calculated  distribution,  this  would  put  too 
much  confidence  in  the  outliers  which  arise  from  using  only  three  points  to  define  the  uncertainty 
in  each  error  source.  Instead  a  scheme  was  developed  to  construct  a  log  normal  cdf  (cumulative 
distribution  function)  from  the  3-point  output  distribution. 

First,  from  the  3-point  output  distribution,  the  median  and  quartilcs  are  determined  (call 
them  d50,  d25,  and  d75).  The  quartilcs  are  used  to  avoid  potential  outliers,  and  the  median  is  used 
instead  of  the  mean  since  the  median  is  a  more  robust  statistic  (small  changes  in  input  values  have  a 
small  effect  on  robust  statistics). 


Secondly,  the  log  normal  parameters — median  (p),  dispersion  (p) — arc  calculated. 


Let  |i 

=  d50 

(3) 

Pu 

=  1  In  (d75/d50)l 

(4) 

Pl 

=  1  In  (d50/d25)l 

(5) 

P 

=  V  pu  •  Pl 

(6) 

The  log  normal  cdf  is  then  defined  so  that  the  5th,  50th,  and  95th  percentiles  can  be  easily 

found. 


F(d)  = 


exp(-x2/2)  dx 


(7) 


where  z(d)  =  ln(d/|i)  /  P 


(8) 
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2.2  DE  ANALYSIS. 


Outputs  from  the  DTE  analysis  include  three  points  describing  the  uncertainty  in  DTE  due 
to  uncertainty  in  EPW  Impact  velocity,  geology,  target  hardness,  tmd  ground  shock  model;  relative 
contributions  of  each  error  source  to  the  uncertainty  in  DTE;  and  a  probability  of  successful 
penetration  (PSP)  which  considers  both  rebounding  and  excessive  axial  loading.  Both  the 
calculated  uncertainty  in  DTE  and  the  calculated  value  for  PSP  are  used  as  inputs  to  the  second  part 
of  the  tarjiet”'"'  uncertainty  methodology — DE  analysis. 

The  analysis  is  concerned  with  the  calculation  of  the  targeting  prediction  measure  DE 
(damage  expectancy),  and  the  evaluation  and  ranking  of  the  error  sources  on  the  uncertainty  in  DE. 

Before  discussing  the  potential  error  sources  treated  in  the  DE  analysis  and  the  details  of  the 
DE  analysis,  some  of  the  major  targeting  prediction  measures  will  be  defined  to  ensure  a  common 
vocabulary. 

Damage  expectancy  (DE)  is  the  probability  of  achieving  at  least  the  specified  level  of  target 
damage,  and  is  the  product  of  two  probabilities — probability  of  arrival  and  probability  of  damage. 

Probability  of  arrival  (PA)  is  the  probability  associated  with  the  weapon  successfully 
arriving  and  detonating  in  the  target  area  as  planned.  PA  accounts  for  pre-launch  survivability  of 
the  delivery  platform,  penetration  of  defenses,  weapon  system  reliability,  successful  penetration, 
and  other  factors  depending  on  the  specific  delivery  platform. 

Probability  of  damage  (PD)  is  the  probability  of  achieving  at  least  the  specified  level  of 
target  damage,  assuming  arrival  and  detonation  of  the  weapon  in  the  target  area.  PD  is  calculated 
using  a  modified  version  of  the  mathematics  in  PDCALC-4  (the  official  targeting  algorithm  used  by 
the  Joint  Strategic  Target  Planning  Staff,  JSTPS).  The  modifications  allow  for  the  explicit 
treatment  of  target  depth,  and  incorporate  the  DUGlC-defined  iso-stress  contour  geometry.  In 
order  to  explicitly  treat  target  depth,  the  distance  damage  prediction  function  (as  shown  in  Figure  7) 
is  defined  in  terms  of  the  slant  range  instead  of  the  commonly  used  ground  range.  The  distance  » 

damage  prediction  function  defines  the  target  damage  probability  when  the  actual  weapon  impact 
location  is  known.  PD  is,  then,  the  integral  over  all  potential  weapon  impact  locations  of  the 
distance  damage  prediction  function. 
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Figure  7.  Calculation  of  probability  of  damage. 


Potential  error  sources  treated  in  the  DE  analysis  include  DTE  (calculated  in  the  DTE 
analysis),  target  depth,  target  location,  target  size,  weapon  aiming  accuracy,  and  distance  damage 
prediction  model.  Each  error  source  is  described  in  the  three-point  format.  Uncertainty  in  target 
location  is  measured  horizontally,  while  uncertainty  in  target  depth  is  measured  vertically. 

The  uncertainty  in  the  distance  damage  prediction  model  is  described  in  terms  of  the  slant 
range  damage  sigma — a  slight  modification  of  the  definition  of  the  damage  sigmo  used  in  the  target 
vulnerability  system.  The  slant  range  damage  sigma  is  a  measure  of  how  rapidly  the  distance 
damage  prediction  function  degrades.  Figure  8  shows  the  effect  of  uncertainty  in  the  slant  range 
damage  sigma  cn  the  distance  damage  prediction  function. 

If  uncertainties  are  assigned  to  each  of  the  six  potential  error  sources,  the  DE  analysis 
would  calculate  729  (3  6)  DE  values — one  for  each  combination  of  error  source  values.  The 
calculated  values  for  DE  are  sorted  to  form  a  cdf.  The  nominal  (planning  value),  minimum, 
maximum,  median  and  quartile  values  for  DE  are  reported,  along  with  one  sigma  type  bounds  for 
the  uncertainty  in  DE. 
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Figure  8.  Effect  of  slant  range  damage  sigma  on  the  distance  damage  prediction  function. 


The  5th  and  95th  percentiles  are  not  of  interest  as  measures  of  the  uncertainty  in  DE  since 
they  are  prone  to  oudiers  and  often  span  the  entire  range  space  for  DE  (from  0  to  1).  Much  more 
informative  are  one  sigma  type  bounds,  the  average  absolute  deviation  about  the  nominal  DE  value, 
calculated  by  the  foilox^ing  equations. 


DEu=  X 

LEi>DEN  ^ 


(9) 


DEl=  X 

DEt<DEM 


DEj  -  DEn 
n* 


(10) 


where  DEu  is  the  average  absolute  deviation  above  the  nominal  DE  value,  DEn  is  the  nominal  DE 
value  found  by  calculating  DE  with  all  error  sources  at  their  nominal  values,  n'  is  the  number  of 
DE  values  greater  than  the  nominal  value,  DEl  is  the  average  absolute  deviation  below  the  nominal 
DE  value,  and  n*  is  the  number  of  DE  values  less  than  the  nominal  value. 
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In  addition  to  all  the  mentioned  statistics  which  quantify  the  uncertainty  in  DE  due  to  the 
uncertainty  in  the  potential  error  sources,  the  same  procedure  used  in  the  DTE  analysis  is  used  to 
calculate  the  relative  contribution  of  each  error  source  to  the  uncertainty  in  DE. 

Comparisons  have  been  made  between  the  distribution  of  DEs  calculated  using  the  3-point 
method  with  that  using  Monte  Caiio  techniques  (Ref.  1,  p  1 19-148).  Continuous  distributions 
(uniform,  normal,  chi-square,  reverse  chi-square)  were  assigned  to  selected  error  sources  and 
propagated  through  the  DE  formulation  using  a  Monte  Carlo  sampling  scheme.  The  resulting 
distributions  of  DE  values  were  compared  to  log  normal  distributions  constructed  from  the  3-point 
output  distributions  (where  the  uncertainty  in  each  eiror  source  is  defined  in  the  3-point  format 
using  the  5th,  50th,  and  95th  percentiles  from  the  assumed  distribution).  The  method  for 
constructing  a  log  normal  distribution  from  an  output  distribution  derived  using  the  3-point  metliod 
was  discussed  earlier  in  this  paper.  The  comparisons  of  distributions  exhibited  very  good 
agreement.  It  was  this  study  that  found  the  best  agreement  with  the  Monte  Carlo  output  to  occur 
when  the  three  points  for  each  error  source  are  defined  as  the  5th,  50ih  and  95th  percentiles. 
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SECTION  3 
EXAMPLE  PROBLEM 


3.1  PROBLEM  DEFINITION. 

The  example  problem  involves  the  application  of  a  low  yield  EPW  against  a  target  sited  in 
a  layered  geology  beneath  a  one  foot  layer  of  concrete.  The  EPW  is  assumed  to  have  a  yield  of 
2U0  Kt,  and  to  impact  the  ground  at  a  velocity  of  800  fps.  It  is  desired  to  assess  the  effectiveness 
of  the  EPW  against  targets  over  a  wide  range  of  depths.  The  target  of  interest  is  a  point  target  of 
moderate  hardness. 


Uncertainties  are  assigned  to  the  following  error  sources. 


Impact  Velocity 

Geology  Layer  Penetrability 

Target  Kandness 

Ground  Shock  Model 

CEP  (Circular  Error  Probable) 

Slant  Range  Damage  Sigma 

Target  Depth 


800  fps  ±  10% 

±  10%  (±5%  for  concrete) 

Nomin^  ±  40% 

X  1.5  (in  stress  at  a  given  depth) 
175,  250  ft 
0.3  ±0.1 

±  20%  over  a  wide  range  of  depths 


All  error  bounds  are  defined  in  terms  of  5th,  50th,  and  95th  percentiles,  except  for 
penetrability  and  CEP,  The  bounds  for  layer  penetrability  represent  a  one  standard  deviation 
spread.  CEP  is  expressed  in  terms  of  two  candidate  values.  The  intent  of  varying  target  depth 
over  a  range  of  values  is  to  show  how  the  uncertainty  in  DE  as  well  as  the  major  contributors  vtiry 
with  different  input  values. 


3.2  EXAMPLE  OUTPUT. 

3.2.1  DTE  Analysis  Output. 

Applying  the  PENDEPTH  code  with  the  mentioned  uncertainty  in  geology  and  impact 
velocity  results  in  a  probability  of  successful  penetration  (PSP)  of  1.0  (i.e.,  no  problem  with  either 
rebounding  or  excessive  axial  loads),  and  nine  depths  of  penetration  ranging  from  3.7  to  6.7 
meters.  Using  DUG  1C,  the  best  estimate  for  depth-to-effect  is  460  meters,  with  5th  and  95th 
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percentiles  of  33..  and  637  meters,  The  coiunlMiiU>mi  of  each  entii  souivp  to  ihr  mueitrtihtv  (o 
DTH  arc  shown  in  I'igurc  d, 


l‘iwure  Heliinve  vontnlnmons  to  I  >  fT  mu  enatmv 


3,2.2  Dl'!  Analysis  Output, 

Best  estimates  ami  umerlainiies  ftir  damaite  espectaius  .ue  nIiowu  m  \  in  ,it  « 
function  o(  target  depth.  Hie  uncertainty  IhiuiuIs  CkmesiHuid  u>  a  •>ne  siaiuhud  desiaiiim  i\)h 
Ixuind;  actually  they  represent  the  average  dts|K>rMttn  diut  Ivlosv  tfie  m«»unal  I  *|.  \aUi»’  tho 
allows  the  ts>unds  to  lx?  asymmetric,  and  ensures  that  the  Umnds  siav  vstiinn  the  li  to  I  inu  i\  ,»i  i.m 
svhich  ni'.  is  defined. 

Notice  ihnt  high  Dt'ls  ate  achieved  lot  taigei  depths  less  ihan  imuis  an. I  itt.si  iii. 
uncertainty  Ixmiuls  arc  widest  when  ihe  lx?sc  estimate  ls*t  I'c  is  tx'isseen  i'  ’  tn.l  o  ' 

Also  shown  III  I'igurc  ID  are  the  relative  CisnittlMiiions  lo  ihe  umeitamts  in  l»|  n«i  to*., 
target  depths.  Notice  that  uncertainty  in  target  lepih  is  not  a  iim)oi  eii\«  somse  toi  shallow  i.uiis  o 
hut  increases  m  imponance  as  the  target  gets  vlee|x?r  .N'so,  ii  can  Ix'  seen  that  un.  eoainis  .m  >  I  f 
Ivcomes  less  a  factor  as  target  depth  inv reuses,  amt  linalh  Oil  um  eitainis  is  a  inahn  >  oniiil>Mi,>« 
lor  all  target  ilepihs.  Ihese  conclusions  regaiding  l>ll'  aitvl  hi  unseiiaints  .md  leloisc 
contributions  arc  sirccific  lo  the  example  pixiblem  inputs 


n  ‘  *  »• 

I  ^44««4^*  >i|rkat 
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SK(TI()N  4 

SUMMARY  OBSKRVATIONS 

The  mciluHlulogy  discussed  in  ihis  paper  quantifies  the  et't'cct  of  systematic  errors  on  the 
tarneiing  prediction  measure  DR  (damage  expectancy).  It  evaluates  and  ranks  the  effects  of 
nitceriainties  for  applications  involving  a  single  buried  target  and  a  single  earth  penetrating  weapon. 

The  mctluxiology  makes  minimal  use  of  assumed  distributions.  Instead  it  dc.scribcs  the 
uiK'oi  tainiy  in  all  but  one  of  the  error  .sources  with  three  points.  All  combinations  of  values  for  the 
error  sources  arc  propagated  through  the  methodology,  creating  a  distribution  of  outcomes. 
Uncertainty  in  the  outcome  is  evaluated  and  the  relative  contributions  of  each  error  source  to  the 
uncertainty  in  outcome  is  ciUculatcd  using  an  :uialysis  of  variance  technique. 

Benefits  of  the  methodology  include  its  simple  design,  its  capability  to  run  on  personal 
computers,  its  fast  nin  time,  the  freedom  allowed  in  specifying  uncertainties  in  error  sources,  and 
its  analysis  of  major  contributors  to  outcome  uncertainties.  The  ranking  of  error  sources  can 
provide  insight  into  deeply  buried  targeting  problems,  possibly  affecting  system  design  and  design 
iradi'offs,  revealing  the  need  for  more  research  into  ground  shock  environments,  or  isolating  the 
effect  of  uncertaindes  associated  with  intelligence  data. 
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SECTION  5 

RUNTIME  PROCEDURES 


This  section  explains  me  requirements  of  the  uncertainty  methodology,  and  the  steps 
required  to  mn  the  code.  The  uncertainties  metlKxiology  consists  of  two  codes  —  DTEUA,  the 
first  part  of  the  methodology  which  calculates  uncertainty  in  depth-to-effect;  and  DUGDE,  the 
second  part  of  the  methodology  which  calculates  uncertainty  in  damage  expectancy. 

DTFUA  (Depth-to-effect  Uncertainty  Analysis)  is  written  in  Microsoft  Quick  Basic,  a 
compiled  language.  DUGDE  is  written  in  Microsoft  Fortran,  version  3.3.  Both  codes  will  run  on 
any  IBM  PC-compatible  machine  with  at  least  512  K  RAM.  A  math  coprocessor  is  not  required, 
but  will  improve  execution  time.  Neither  a  hard  disk  nor  graphics  capability  are  required.  DTEUA 
does  require  a  printer. 

Conventions  used  in  these  codes  and  in  this  documentation  will  now  be  discussed.  All 
input  to  be  entered  by  the  user  is  shown  underlined.  Each  command  must  be  followed  by  a  return. 
Options  are  often  given  inside  parentheses  and  separated  by  the  symbol For  example,  "(Y/N)" 
means  that  the  user  should  input  either  a  Y  or  an  H-  Options  are  not  case  sensitive,  i.e.,  the  input 
"Y"  and  "y"  are  Heated  identically.  If  the  user  do^-,  not  enter  a  satisfactory  input,  the  program  will 
prompt  the  user  with  the  same  question.  All  values  entered  by  the  user  are  checked  to  ensure  their 
validity. 


Default  values  are  given  inside  angle  brackets.  For  example,  "(Y/<N>)".  In  this  case, 
"N"  is  the  default  value.  If  the  user  wants  to  input  "Y",  he  must  enter  either  X.  or  y.  Otherwise  the 
default  value  is  used.  The  easiest  way  to  get  the  default  is  just  to  hit  return. 


5.1  RUNNING  DTEUA  (DEPTH-TO-EFFECT  UNCERTAINTY  ANALYSIS). 

Steps  required  to  run  DTEUA  are  listed  below.  DTEUA  does  require  a  printer.  Before 
running  the  code,  the  user  should  make  sure  that  the  DTEUA  disk  or  subdirectory  contains  the 
following  files  required  by  the  program;  DTEUA.EXE,  DUG1C.EXE,  DUGIC.MSG, 
PROPERTY.DAT,  and  a  DUG  1C  input  file  which  specifies  the  geology. 
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Insert  the  DTEUA  disk  into  drive  A:, 

At  the  DOS  prompt,  type  Al  to  set  the  default  drive  to  A:. 

Type  DTEUA  to  begin  execution. 

The  program  will  load  and  will  check  to  see  if  the  printer  is  turned  on  and  is  on-line.  If 
not,  the  user  will  see  a  "testing  printer"  message  on  the  screen.  That  message  will  not 
disappear  until  the  printer  is  re^y. 

A  banner  screen  will  appear  on  the  monitor;  to  continue  the  user  should  hit  the  return  key. 

DTEUA  will  prompt  for  the  user  to  enter  the  filename  to  which  results  should  be  written. 

The  user  will  be  prompted  to  enter  the  DUG  geology  filename.  This  is  one  of  the  files 
read  by  DUG  1C  and  which  must  exist  prior  to  running  DTEUA. 

The  code  prompts 

Save  DUG  output  files?  (<Y>/N): 

If  the  user  wants  to  save  the  files  created  by  DUG  1C,  he  should  hit  return;  otherwise  type 

a. 

The  user  is  asked  for  information  concerning  impact  velocity.  First  the  user  is  asked  if 
‘  uncertainty  bounds  are  desired.  If  so,  the  user  is  asked  to  enter  a  nominal  value,  a  lower 
bound,  and  an  upper  bound  for  impact  velocity.  If  no  uncertainty  bounds  are  desired,  the 
user  is  asked  to  enter  a  single  value  for  impact  velocity.  Each  parameter  has  a  default 
value. 

Similar  information  for  target  hardness  is  requested  next,  along  with  a  single  value  for 
weapon  yield. 

Other  weapon  inputs,  those  required  to  calculate  depth  of  penetration,  aie  prompted  for 
—  including  impact  angle,  penetrator  weight,  penetrator  diameter,  nose  shape  (cone  or 
ogive),  nose  performai.ee  coefficient,  and  nose  length. 

Weapon  inputs  are  followed  by  geology  inputs.  After  the  user  specifies  the  number  of 
layers  (between  1  and  15),  a  worksheet  for  entering  geological  information  by  layer 
appears  on  the  screen.  Thickness  in  feet,  S-number,  and  layer  type  (soil,  rock,  concrete, 
or  fi-ozen  soil)  must  be  entered  for  each  layer. 

ITie  user  may  use  the  TAB  key  or  any  of  the  keyboard  arrows  to  move  the  highlighted 
cell  about  the  screen.  All  inputs  goes  into  this  cell.  After  entering  a  value,  move  to  the 
next  cell  by  pressing  return,  TAB,  or  an  arrow  key. 

To  change  tlie  layer  type,  position  the  highlighted  cell  anywhere  on  the  desired  layer,  then 
press  the  space  bar  to  toggle  to  layer  type  from  soil  to  rock  to  concrete  to  frozen  soil  back 
to  soil.  As  the  user  toggles  the  layer  type,  default  values  for  that  layer's  S-number  may 
appear  in  the  correct  cell  on  the  screen.  The  default  value  for  concrete  is  1,  and  for  frozen 
soil  it  is  1.7.  A  note  about  the  layer  types;  As  mentioned  previously,  there  is  a  constraint 


as  to  which  layers  may  be  frozen.  Normally  only  the  top  layer  may  be  frozen.  The 
exception  is  when  the  top  layer  is  concrete;  only  then  can  the  second  layer  may  be  frozen. 
If  the  layer  cannot  be  frozen,  the  space  bar  will  toggle  the  layer  type  from  soil  to  rock  to 
concrete  back  to  soil,  skipping  the  frozen  soil  layer  type. 

The  user  may  not  place  the  highlighted  cell  at  the  thickness  for  the  bottom  layer.  It  is 
assumed  that  this  tliickness  is  infinite.  If  in  doubt  about  an  S-number,  the  user  may  type 
ll  (note  the  lower  case)  for  a  help  screen.  After  typing  ll,  the  user  may  hit  any  key  to 
return  to  the  input  screen.  If  the  user  types  an  invalid  key,  the  program  beeps. 

After  entering  thicknesses  and  S-numbers  for  all  layers,  the  user  should  hit  the  escape 
key  (ESC).  All  values  are  checked  to  make  sure  they  are  valid.  If  any  errors  arc  found, 
the  input  screen  reappears,  complete  with  error  messages,  and  with  those  cells  in  error 
highlighted. 

In  addition  to  nominal  values  for  layer  thickness  and  penetrability  (S-number), 
uncertainty  values  are  also  entered.  First  the  user  is  asked  for  the  number  of  desired 
Monte  Carlo  trials  —  at  least  20  must  be  specified  so  that  the  5th  percentile  may  be 
found. 

Then  a  table  of  nominal  thicknesses  and  S-numbers  is  printed  to  the  screen.  The  user 
must  enter  lower  and  upper  bounds  for  thickness  and  a  one  standard  deviation  for  the  S- 
number  uncertainty.  If  no  uncertainty  is  desired  for  a  given  layer  thickness  or  S-number, 
die  user  should  hit  the  return  key.  If  no  uncertainty  is  chosen  for  any  layer  (thickness 
and  S-number),  the  code  will  perform  only  one  Monte  Carlo  trial  regardless  of  the 
number  entered  by  the  user. 

Lastly,  the  remaining  inputs  pertain  to  calculating  PSP  (probability  of  successful 
penetration)  —  includuig  constraints  on  A-ebounding  depdi  and  maximum  axial  loading. 

The  program  will  ensure  that  all  necessary  files  are  available;  if  not,  it  will  notify  the  user 
of  the  missing  file.  The  designated  number  of  Monte  Carlo  trials  will  be  run  for  each 
specified  impact  velocity.  DUG  1C  will  then  be  run  three  times  for  each  impact  velocity  , 
once  for  each  if  no  uncertainty  was  assigned  to  the  geological  information,  After 
DUG  1C  is  run,  its  output  is  examined  to  determine  depths-to-effect.  Relative 
contributions  of  each  error  soiuoc  are  calculated,  along  with  a  PSP  value,  and  predicted 
uncertainty  in  DTE.  An  output  file  (previously  named  by  the  user)  is  created  and 
execution  stops. 


5.2  RUNNING  DUGDE  (DAMAGE  EXPECTANCY  UNCERTAINTY 
ANALYSIS). 


Required  steps  to  run  the  second  part  of  the  uncenainties  methodology,  DUGDE,  are 
explained  below.  No  files  other  than  the  executable  (DUGDE.EXE)  are  required.  DUGDE  may 
be  run  after  running  DTEUA;  or  if  a  DTE  is  provided  by  some  other  means,  DUGDE  may  be  run 
without  any  input  piovided  by  DTEUA. 
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Insert  the  DUGDE  disk  into  drive  A:. 


•  At  the  DOS  prompt,  type  Ai  to  set  the  default  drive  to  A:. 

•  Type  DUGDE  to  begin  execution. 

•  The  program  will  load  ?r  d  will  prompt  the  user  if  output  should  go  to  the  printer  or  to  the 
screen.  Tire  user  should  enter  £  for  printer,  the  screen  is  the  default. 

•  The  next  input  is  if  the  target  is  a  point  or  an  area  target.  Area  targets  are  assumed  to  be 
circular  and  to  follow  a  circular  normal  distribution,  where  size  is  measured  in  terms  of 
an  R95  (the  radius  of  a  circle  which  encompasses  95%  of  the  target  elements). 

’  Next  the  user  is  asked  whether  the  weapon  detonates  at  deptli  of  burst  (DOB).  If  so,  the 

user  is  asked  to  supply  a  DOB  in  feet.  The  user  is  not  asked  to  specify  an  uncertainty  in 
DOB  —  this  effect  is  taken  into  account  by  the  depth-to-effect  parameter.  The  DOB  is 
used  to  calculate  an  "effective"  target  depth,  the  vertical  distance  between  weapon  and 
target  equal  to  the  actual  target  depth  minus  the  DOB.  In  addition,  when  the  user 
specifies  a  DOB,  a  different  iso-stress  contour  is  used  compared  to  that  used  for  a  surface 
burst  The  difference  in  the  contours  is  illustrated  in  Figure  1 1. 


Contact  Burst 


Depth  of  Burst 


Figure  11.  Two  DUG  1C  iso- stress  contours. 


'fhe  next  four  inputs  pertain  to  the  probability  of  arrival.  ITie  user  is  asked  to  enter  values 
for  PLS  (pre-launch  survivability),  FTP  (probability  to  penetrate),  WSR  (weapon  system 
reliability),  and  PSP  (probability  of  successful  penetration).  Each  is  a  factor  between  0 
and  1.  PSP  is  calculated  by  DTEUA. 

The  user  is  asked  to  hit  return  to  continue. 

The  3-point  parameter  input  screen  appears.  For  more  information  on  the  3-points,  see 
Section  2. 
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The  user  may  use  the  TAB  key  or  any  of  the  keyboard  arrows  to  move  the  highlighted 
cell  about  the  screen.  All  input  goes  into  this  ceU.  As  the  user  moves  the  cursor  around 
the  screen,  the  program  prints  a  three  line  message  describing  the  parameters  at  the 
bottom  of  the  screen. 

Nominal  values  must  be  entered  for  all  parameters.  There  is  no  requirement  that  both 
upper  and  lower  values  be  input,  or  that  bounds  be  symmetric  about  the  nominal  value, 
^en  the  user  is  satisfied  with  the  inputs,  he  should  move  the  highlighted  cell  down  to 
the  bottom  of  the  worksheet,  next  to  the  "Accept?:"  prompt  (alternatively  hit  the  end  key). 
Then  type  either  Y  to  continue  or  M  to  reinitiate  the  input  process.  The  return  key  should 
not  be  hit.  If  the  user  enters  R  he  will  be  instructed  to  enter  all  inputs  again.  Otherwise, 
the  program  will  check  to  make  sure  that  all  parameters  have  nominal  bounds,  and  that  all 
inputted  lower/upper  bounds  are  indeed  less  than/greater  than  the  nominal  value.  Any 
inconsistencies  or  missing  values  will  be  brought  to  the  user's  attention  by  flashing 
parameter  names.  The  user  must  correct  his  input  and  again  move  the  highlighted  cell  to 
the  "Accept?:"  prompt  and  enter  Y- 

The  program  will  perform  its  calculation  and  print  the  output  to  the  screen  or  to  the 
printer,  as  specified  by  the  user.  If  the  output  device  is  the  screen,  the  user  may  have  to 
hit  the  return  key  to  see  the  Figures  of  Merit  chart.  (The  Figures  of  Merit  chart  is  only 
printed  if  at  least  one  bounding  value  was  entered.) 

The  program  prompts  whether  the  user  would  like  to  write  the  data  to  disk.  If  the  user 
enters  Y.  he  will  then  be  asked  to  enter  the  save  file  name. 

The  program  prompts  the  user  for  more  calculations.  If  additional  calculations  are  not 
required,  the  user  should  type  M*  Otherwise,  type  Y;  execution  will  begin  again,  and  all 
previously  inputted  3-point  parameter  values  wUl  be  remembered. 


I 
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